Universal transport in 2D granular superconductors 
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The transport properties of quench condensed granular superconductors are presented and ana- 
lyzed. These systems exhibit transitions from insulating to superconducting behavior as a function 
of inter-grain spacing. Superconductivity is characterized by broad transitions in which the resis- 
tance drops exponentially with reducing temperature. The slope of the log 

R versus T curves turns out to be universaly dependent on the normal state film resistance for all 
measured granular systems. It does not depend on the material, critical temperature, geometry, or 
experimental set-up. We discuss possible physical scenarios to explain these findings. 
PACS: 74.80.Bj, 74.76.Db, 74.50.+r 
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Two dimensional granular superconductors, i.e. sys- 
tems of superconducting grains embedded in an insulat- 
ing matrix, exhibit a superconductor to insulator transi- 
tion (SIT) as a function of the inter-grain distance. Such 
systems are of interest since they can serve as a model 
system for many dirty superconductors. In particular it 
has been suggested that a granular superconductor can 
mimic the behavior of high Tc superconductors and its 
properties can be used to explain some of the features 
observed in the cuprates 

An established technique to study the properties of 
granular superconductors is quench condensation 
In this method one performs sequential evaporations on 
a cryogenically cold substrate under UHV conditions us- 
ing the following scheme: Metallic leads for four-terminal 
measurements are prepared on an insulating substrate 
which is then mounted onto an evacuated measurement 
probe and immersed in a liquid He bath. The low tem- 
perature of the probe causes cryopumping and hence the 
substrate is situated in UHV conditions and, at the same 
time, at cryogenic temperatures. This enables evapora- 
tion of ultra-clean thin superconducting layers on a sub- 
strate held at temperatures lower than 10K while con- 
tinuously monitoring the film resistance and thickness. 
Once a desired resistance (or film thickness) is achieved, 
the evaporation is terminated and the transport prop- 
erties are measured. Incremental layers of material are 
then added in-situ and further measurements are taken 
at different film resistances. Using this method one can 
study the properties of a single sample throughout the 
entire transition from an insulator to a superconductor 
as a function of the amount of deposited material while 
keeping the sample at low temperatures and in a UHV 
environment without having to thermally cycle it (risking 
metallurgical or structural changes due to annealing) or 
to expose it to atmosphere (thus oxidizing the surface). 

If the samples are quench condensed on a passivated 
substrate (such as SiC^), they grow in a granular manner 
so that the film begins its growth as disconnected islands 
with diameters of 100-200 A |g-|]. The average distance 
between the islands decreases upon adding material. Be- 



yond a percolation threshold, the grains connect, forming 
a continuous conducting layer. In these samples there is 
a critical nominal thickness, dc, below which no con- 
ductivity can be measured (the sheet resistance, R, is 
larger than 10 10 S1). Once the thickness, d, of the sample 
is larger than dc, R drops exponentially with thickness 
until, for R < 6k£l, it crosses over to a normal ohmic 
behavior (Ral/d) @. 

Varying the thickness of the film causes a transition 
from an insulating behavior for the thinnest films to 
a superconducting behavior for thick films. Figure 1 
demonstrates examples for this transition in three differ- 
ent granular superconductors: a Pb film (critical temper- 
ature, T c , -7.2 K), a Sn film (T c «4.5K) and a Pb/Ag 
bilayer. The latter is a system in which a thin layer 
of insulating granular Pb is quench condensed on a Si02 
substrate followed by sequential evaporations of Ag ultra- 
thin layers . The curves were measured in a shielded 
room using standard 4 probe technics and assuring, for 
each point on the curve, that the I-V characteristics are 
in the linear regime. 

The three systems show a transition from an insulator 
for thin films to a superconductor at thicker films. In all 
2D granular systems these transitions are characterized 
by broad resistance tails and critical temperatures which 
are not very well defined until the resistance of the film 
becomes low. The transitions become sharper as mate- 
rial is added to the film. In this paper we define Tc as 
the temperature at which the resistance starts dropping 
exponentially with lowering the temperature. We justify 
this by observing that this is the temperature in which 
the individual grains become superconducting even in the 
insulating case as discussed below. In the Pb and Sn films 
(top two graphs in figure 1) Tc has bulk value and barely 
changes throughout the entire transition. Moreover, even 
on the insulating side of the SIT the curve changes its 
slope at T = Tc reflecting the presence of superconduc- 
tivity even in the thinnest measurable samples. The fact 
that the grains are superconducting with bulk properties 
throughout the entire transition has been demonstrated 
by tunneling measurements |10| . A bulk energy gap was 
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observed in the grains even when the film was on the 
insulating side of the SIT. 
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FIG. 1. Resistance versus temperature for sequential lay- 
ers of quench-condensed granular Pb (Top), Sn (middle) and 
Pb/Ag (bottom). Different curves correspond to different 
nominal thickness. 



In the Pb/Ag film the Tc decreases as material is 
added to the film. This is due to suppression of supercon- 
ductivity in the Pb grains by the silver overlayer because 
of proximity effects. Nevertheless, the general behavior of 
the insulator to superconductor transition is very similar 
to that of the pure superconductors. 

The transport curves shown in figure 1 are attributed 
to the unique nature of the SIT in a granular structure 
pr| . In these systems each grain sustains superconduc- 
tivity with bulk properties. However, for the high resis- 
tance samples, there are rapid, thermally activated phase 
fluctuations between the grains leading to an insulat- 
ing R-T curve. As the grains become closer the inter- 
grain Josephson coupling increases until phase coherence 
is achieved throughout the film. 

A striking feature common to the quench condensed 
granular superconductors is that at temperatures be- 
low Tc the sheet resistance decreases exponentially with 
decreasing temperature and can be described approxi- 
mately by an "inverse arrennius" law: 

R = R e^ (1) 

It is important to note that we see no flattening of the 
R-T curve ( resistance approaching a flat temperature 
dependence at low temperatures) in any of our samples 
and in no regime of thicknesses |12]] . Nor does the R-T 
change its trend and begin dropping to zero resistance at 
low temperature. The dependence described in equation 
1 is observed in all of our granular samples using differ- 
ent materials and different measurement apparatuses. It 
spans many orders of magnitude in R, persists to tem- 
peratures below 100 mK |ll|] and is observed for the dif- 
ferent steps of the sequential evaporations providing the 
resistance is smaller than a few tens of fcfi. Moreover, 
the slope of the log R versus T curves, 1/T of equation 
1, turns out to be universal for all of our samples. It 
depends only on the normal state sheet resistance, Rn, 
and does not depend on the material, the critical temper- 
ature or sample geometry. We illustrate this in Figure 2 
where we show the dependence of the inverse slope, To, 
on Rn for a large number of granular superconductors 
(different materials, geometries and prepared in various 
quench condensation evaporators) . It is seen that all the 
slopes fall on a master-plot having the form of: 

T « C*R N (2) 

where C* is a constant of approximately 0.05-^-. Note 
that we have included samples having different critical 
temperatures as well as the Pb/Ag system in which Tc 
varies with the thickness (and hence with Rn). 

This observation indicates that the behavior of the su- 
perconducting tails does not depend on the properties 
of the superconducting grains themselves but only on 
their geometrical arrangement (density, inter-grain spac- 
ing, morphology configuration etc.) which determines 
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the tunneling percolation network and Rn- Two super- 
conductors having the same Rn but different critical tem- 
peratures Tci and Tqi will have parallel superconducting 
R-T tails, at temperatures lower than the respective Tc, 
shifted by Tqi — Tc2- An example for such behavior is 
shown in the insert of figure 2. 



are expected to be Josephson coupled with a Josephson 
binding energy of the form: 




FIG. 2. The inverse slope of the R-T curves ( To of equa- 
tion 1) as a function of the normal-state sheet resistance for 
different granular samples. The open symbols are the ex- 
perimental data. Squares are results for various Pb sam- 
ples, circles are for Sn samples and triangles for Pb/Ag 
proximity samples. Single symbols are used for the same 
quench-condensed films at different thicknesses (and hence 
Rjvs). Full symbols are the slopes extracted from the calcu- 
lated samples of figure 3 (squares -Pb and circles - Sn). The 
line is a guide to the eye. Insert: Resistance versus tem- 
perature for granular Pb (squares), Sn (circles) and Pb/Ag 
(triangles) samples having a normal state resistance of 4 kfi. 

We can now combine this observation with equations 
1 and 2 and extract the expression for Rq of equation 1: 



Ro = Rnz t o = R N e~' 



(3) 



Ro is the value of the resistance obtained by extrap- 
olating the R-T curves to T=0. This means that if the 
observed exponential R-T curves persist to zero temper- 
ature the samples will exhibit a finite zero-temperature 
resistance which will depend both on R^ and on Te- 
la an attempt to examine possible intuitive under- 
standing of the observed phenomena we note that at tem- 
peratures below bulk Tq the grains have been shown to 
be fully superconducting |1C( ]. Hence, each two grains 
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where A(T) is the temperature dependent supercon- 
ducting gap and Rat is the normal resistance between 
the grains. One can expect that the ratio ^fi T would 
determine whether the grains are phase coupled or not. 
As the temperature of a granular system is lowered Ej be- 
comes larger than K^T in an increasing number of pairs 
of grains thus increasing the effective superconducting 
regions which are phase coherent. This leads to a char- 
acteristic "phase-coupled" length which grows with de- 
creasing temperature until superconducting percolation 
is achieved. 

However, the naive model described above can not ex- 
plain the observed broad superconducting tails in our 
films. Since the samples are two dimensional, scaling 
considerations imply that as long as superconducting per- 
colation is not achieved throughout the sample, the re- 
sistance of the granular system should be independent 
on the size of the superconducting clusters. The growth 
of the phase coherent clusters with decreasing tempera- 
ture would not change the overall 2D resistance. One can 
expect to see a flat resistance versus temperature curve 
below the percolation threshold and a sharp supercon- 
ducting transition at the percolation threshold. The fact 
that the experimental transitions are broad implies that 
the resistance of the individual junctions is temperature 
dependence. This urges us to consider superconductor 
fluctuation effects as fluctuations introduce a broad tem- 
perature dependence in a Josephson junction resistance 
rather than a sharp transition at Ej=KsT. 

The temperature dependence of the zero-bias resis- 
tance of each pair of grains due to thermal phase fluc- 
tuations is expected to take the Ambegaokar-Halperin 
form liiJlSl: 



R{T) = 



Rn 



('o(W)) 2 



(•5) 



where Io(X) is the modified bessel function of order 
0. Using this model we have calculated the temperature 
dependence of ID arrays of junctions subject to thermal 
fluctuations. The results of the simulations for Pb and 
Sn samples having different Rn are shown in figure 3. It 
is seen that R(T) may qualitatively mimic an exponen- 
tial dependence for temperatures larger than about 
As Rn is reduced the R-T curves become sharper in a 
similar way to that seen in the experiments. While the 
curves are not truly e T dependent, for purposes of illus- 
tration we can extract an approximated slope from each 
curve. The extracted slopes are plotted on the master- 
plot of figure 2 showing relatively good agreement with 



3 



the experimental results. We can not expect better than 
qualitative agreement as this is a simple ID chain model. 
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FIG. 3. Numeric simulations of the resistance of arrays of 
10 junctions in series based on equation 5. The top graph is 
for Sn and the bottom for Pb. 

Hence, the thcrmal-fluctuation-picture can provide 
some insight to the physics of the 2D granular films. 
However, two major points in which the simple model 
deviates from the experimental data have to be noted: 
In the first place the theoretical model predicts that 
the slopes should depend on Our simulations fail 
to reproduce the Tc independence and the "universal- 
ity" of the experimental transition tails. The calculated 
curves for Sn and Pb samples having the same Rn yield 
Sn slopes which are consistently sharper than those of 
Pb. Secondly, the thermal-fluctuation model can qual- 
itatively imitate the experimental behavior only for a 
limited range of temperatures. At low temperatures the 
thermal fluctuations are expected to freeze out causing 
the resistance to drop much more rapidly, until, for zero 
temperature, the resistance reaches zero ( Ro=0). We do 
not see signs for such a tendency in the experiments even 
at very low temperatures. We have considered the possi- 



bility that at low temperatures quantum phase fluctua- 
tions dominate the behavior and induce a finite resistance 
even at T=0. The presence of such effects is consistent 
with the observed experimental dependence of Rq on Tc 
as shown in equation 3 since the quantum fluctuations 
are expected to decay exponentially with Ej. However, 
we are not able to reproduce the experimental exponen- 
tial R-T curves over the entire temperature range using 
a model which combines thermal fluctuations and quan- 
tum fluctuations. Such a model would suggest that at 
low temperature the R-T curve would flatten out and 
saturate at a constant value. As noted above, we do not 
observe such behavior in any of our samples. Clearly, 
further theoretical treatment is required in order to shed 
more light on the origin of the findings presented in this 
paper. 
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